THE STRATONOVICH HEAT EQUATION: A CONTINUITY RESULT 
AND WEAK APPROXIMATIONS 
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Abstract. We consider a Stratonovich heat equation in (0, 1) with a nonhnear multipHcative 
noise driven by a trace-class Wiener process. First, the equation is shown to have a unique 
mild solution. Secondly, convolutional rough paths techniques are used to provide an almost 
sure continuity result for the solution with respect to the solution of the 'smooth' equation 
obtained by replacing the noise with an absolutely continuous process. This continuity result 
is then exploited to prove weak convergence results based on Donsker and Kac-Stroock type 
approximations of the noise. 



1. Introduction and main results 

The main motivation of the paper comes from |3j, where the authors consider, for some fixed 
T > 0, the stochastic heat equation 

^(*'^)-^(*'^) = ^"(^'^)' (t,^)G[o,r]x[o,i], (1) 

with some initial data and Dirichlet boundary conditions, where the random fields (0")n>i 
verify that the family of processes 6^{t,x) := 6^{s,y) dyds converge in law, in the space 
C([0,T] X [0, 1]) of continuous functions, to the Brownian sheet. Then, sufficient conditions on 
6"' are provided such that converges in law, as n — )• oo, to the mild solution Y of 

dY d'^Y 

where W{t, x) stands for the space-time white noise. Applications of this result include the 
case of a Donsker type approximation, as well as a Kac-Stroock type approximation in the 
plane. 

Such diffusion approximation issues for stochastic PDEs have been extensively studied in the 
literature. Let us quote here Walsh [32], Manthey [2D1I2I], Tindel [3U], Carmona and Fouque 
[S], Florit and Nualart [T2], just to mention but a few. 

Now, following the line of [3], a natural question to be dealt with is to try to get the same 
type of weak convergence in a non-additive situation, that is when the term 9^{t,x) in ([1]) is 
replaced with f{Y^{t,x))9'^{t,x), for some sufficiently smooth function / : M — >■ M. In this 
case, one expects that the limit equation is of Stratonovich type, as it was the case in [8] and 
|12j (see also [21 [29] for examples of a similar behaviour). This phenomenon has been recently 
illustrated by Bal in [1] as well, for the weak approximation of a linear parabolic equation in 
with random potential given by Y^{t, x)d"'{x). 
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Going back to our setting, and focusing first on what we expect to be our limit equation, 
we should consider the Stratonovich heat equation 

^{t,x)-^{t,x) = f{Y{t,x))oW{t,x), (t,x) G [0,r] X [0,1]. (2) 

Unfortunately, a well-known drawback in this situation is that the solution admits only very 
low regularity (see |31j). a major obstacle for our treatment of the non-linearity of the problem. 
For this reason, we have chosen to restrict our attention to the case of a trace-class noise. To 
be more specific, we will assume that W is the formal derivative of a L^(0, l)-valued Wiener 
process {Wt, t G [0, T]} with covariance operator Q satisfying the following property: 

Hypothesis 1. Let {ek)k>i be the basis of eigenf unctions for the Dirichlet Laplacian A in 
1/^(0,1) given by ek{x) := ^/2sul{k^Tx), x G [0,1]. We assume that there exists a sequence 
of non-negative real numbers {Xk)k>i ^^'d a parameter rj > such that Qe^ = ^k^k for every 
k > 1 and '^k>ii^kk^^) < oo. Without loss of generality, we assume that rj G (0, |). 

In particular, for any fixed t > 0, the process Wt can be expanded in L^(il; L^(0, 1)) as 

Wt = Y,V^k/3^^k, (3) 

k>l 

where {I3^)k>i is a family of independent Brownian motions. Note that the condition X]fc>i(Afc • 
k^'^) < oo is only slightly stronger than the usual trace-class hypothesis Ylik>i < oo, insofar 
as r] can be chosen as small as one wishes. For instance, it covers the case where Q = (Id — A)~^ 
with r > ^. 

Another change with respect to [3] lies in our formulation of the study: compared to the 
random field approach in [3], here it has turned out to be more convenient to use the Hilbert- 
space- valued setting of Da Prato and Zabczyk [9]. In particular, we are interested in the mild 
form of equation ([2]) , which is given by 

Yt = StiP + /* St-u{f{Yu) o dWu), t G [0, T], (4) 
Jo 

where from now on, we will use the notation Yt{-) := Y(t, •), ■0 is some initial condition and 
/ : M — )• M is a smooth enough mapping. As usual, (St)t>o denotes the strongly continuous 
semigroup of operators generated by —A. 

A first part of the paper (Section [2D will be devoted to the interpretation of (jH) as a 
Stratonovich equation, and it will allow us to exhibit an existence and uniqueness result for 
the solution. We should mention here that the stochastic heat equation in the Stratonovich 
framework has already been studied in various settings, most of them in the case of a linear 
multiplicative noise (see e.g. [U 13 |T7]). Once we have given a full sense to ([3]), our strategy 
to study weak approximations of the solution could be stated in the following loose form: 

(a) We will first establish an almost sure continuity result (in some suitable space-time topol- 
ogy) for the solution of (jH) with respect to the solution of the 'smooth' heat equation obtained 
by replacing W with an absolutely continuous process W (see Theorem II. ip . 

(b) Then, for two particular families of absolutely continuous processes approximating W, we 
will rely on our continuity result to show convergence towards the solution in some possibly 
different probability space, leading us to the expected weak convergence (see Theorem ll.2|) . 

Our strategy to compare the solution y of with 'smooth' solutions is based on a genuine 
rough-paths type expansion of the equation, which follows the ideas of [16l HU [10]. Rough- 
paths techniques have indeed proved to be very efficient as far as approximation of non-linear 
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systems in finite dimension is concerned (see \1'6\ Chapter 17]), and it is therefore natural to 
address the same question in this infinite-dimensional background. Note that the model given 
by (HD differs from those studied in [lH llUj. where only finite-dimensional noises are considered, 
forcing us to revise most of the technical details behind this procedure (see Section E]). 

In order to state the above-mentioned results with more details, we need to introduce the 
spaces in which our random variables will take their values. First, as far as the spatial regularity 
is concerned, the fractional Sobolev spaces must come into the picture. Thus, for every a € M 
and p > 2, we will denote by Ba,p the fractional Sobolev space of order a based on L'p{0, 1), 
that is 

¥PGi3„,p^(-A)>GLP(0,l), 
where A stands for the Dirichlet Laplacian in L^(0, 1) (see e.g. [26] for a thorough study of 
these spaces). For the sake of conciseness, we will write Ba for Ba,2 and B for Bq = L'^{0, 1) 
throughout the paper. We will also denote by Boo the set of continuous functions on [0,1], 
endowed with the supremum norm. 

Of course, we will also have to deal with the time regularity of our processes. So, for any 
subinterval / C [0,T] and any Banach space V, we define C^{I; V) as the space of continuous 
functions y : I ^ V and set 

M[y;C^iI;V)] ■.= sup\\yt\\v. 

t&i 

Moreover, for any A > 0, we introduce the space C"^ (/;!/) of A-Holder continuous valued 
functions endowed with the seminorm 

M[y;C\l;V)]:= sup M^Ifjk. (5) 

s<tei \t — s\ 

Note that in the case where / = [0, T], we will often write C^{V) for C'^([0,T]; V). 

Now, consider any process W defined on the same probability space as W and with ab- 
solutely continuous paths in Bri^2p, for every integer p > I (recall that rj has been defined in 
Hypothesis [1|) . Then, let {If, t £ [0,7"]} be the unique solution of the Riemann-Lebesgue 
equation (considered in a pathwise sense): 

Yt = St^ + /* St-u{f{Yu) ■ dWu), t e [0, T], (6) 
Jo 

where tp £ B. As evoked earlier, our first main result will consist in comparing such a solution 
Y with the solution y of Q. This result can be stated as follows. 

Theorem 1.1. Assume that Hypothesis [J\ holds true for W and some parameter r] > 0, and 
let f : W ^ be a function of class C^, bounded with bounded derivatives. In addition, pick 
7 G (^, ^ -|- 7?) and assume that the initial condition V' (resp. il^) in ^ (resp. in ^) belongs 
to Bry . Then there exist e > and p > 1 such that 

M[Y-Y;C\B^)]<F,,p(^mB,AmB,M[W;d-'{Br,,2p^^^ 

{ U - + M[W - W- C^^(e,,2p)] } , (7) 

for some deterministic function F^^p : (]R+)^ — t- M+ bounded on bounded sets. 

The topologies involved in this statement are directly inherited from our rough-paths analysis 
of the equation, and their relevance should therefore become clear through the lines of Section 
[3] (see in particular the proof of the central Proposition 13. 9p . Note that this bound certainly 
remains valid with respect to some Holder norm (in time) for the left-hand side of ([7]), as our 



4 



AURELIEN DEYA, MARIA JOLIS AND LLUIS QUER-SARDANYONS 



arguments will suggest it. However, due to the technicality of the rough-paths procedure, we 
have preferred to focus on the behaviour of the supremum norm (see also Remark l3.12|) . 

Our next step will consist in applying the above Theorem ll.il - on some possibly larger prob- 
ability space - to two particular families of absolutely continuous processes that approximate 
W, so as to retrieve weak convergence results for the solution. To define these approximation 
processes, we will make use of the following additional notation. Namely, on a probability 
space {Q,J-,P), given a sequence {X^)k>i of centered i.i.d processes admitting moments of 
any order, we set, for every t > 0, 



Thanks to our forthcoming Proposition 14.31 we know that W(X ) is indeed well-defined as 
a process on {Q,J^,P) with values in 13rj^2p, for all p > 1. Let us also specify that, given a 
sequence (/3")n>i of real- valued processes, we will henceforth denote by (/?""'%>! = {P"''^)n,k>i 
a generic sequence of independent copies of (/3"')n>i (defined on a possibly larger probability 
space) . 

The two families of approximations at the core of our study can now be introduced as follows 
(we fix T = 1 for the sake of clarity): 

(i) The Donsker approximation W"" := W(5'^''), where is a sequence of appropriately 
rescaled random walks. To be more specific, let {Zj)j>i be a family of i.i.d random variables 
with mean zero, unit variance and admitting moments of any order. Then, for each n G N, set 



Recall that, by Donsker Invariance Principle (see e.g. |19l Thm. 4.20]), S"" is known to converge 
in law to the standard Brownian motion in C'^([0, 1];M), as n — t- oo. 

(ii) The Kac-Stroock approximation := W(0"' ), where 0" stands for the classical Kac- 
Stroock approximation of the (one-dimensional) Brownian motion. Precisely, introduce a stan- 
dard Poisson process N and a Bernoulli variable C independent of A^, with P[<^ = 1) = 1/2. 
Then, set 



Here again, the sequence (6'")„>i thus defined converges in law in ^^([0, 1];M), as n — )• oo, to 
a standard Brownian motion (see e.g [18\ I24j). 

Of course, the one-dimensional weak convergence of 5*^ (resp. 0") towards the Brown- 
ian motion is a priori not sufficient for us to apply Theorem 11.11 Our aim is to turn this 
one-dimensional weak convergence into an almost sure convergence result for W(5'"'') (resp. 
W(0"'')) with respect to the topology involved in ([Tj), and this will appeal in particular to 
Skorokhod embedding arguments (see Section 2]) . Together with Theorem 11.11 the strategy 
ends up with the following statement. 

Theorem 1.2. Under the hypotheses of Theorem \ fix an initial condition ip = tp ^ B^, 
and denote by the (Riemann-Lehesgue) solution of ^ associated with either the Donsker 
approximation = W(5"''') or the Kac-Stroock approximation = W{9"'''). Then, as 
n — 7- oo, converges in law to Y in the space C^(B-y). 

The paper is organized as follows. Section [2] is devoted to a few preliminaries on the theoreti- 
cal study of the Stratonovich heat equation ([2]). The rough-paths type analysis of this equation 



W{X-)t ■.= Y,^kX^ek. 



k>l 





(9) 
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is performed in Section O and it will lead us to the proof of our continuity result Theorem 11.11 
In Section HI we will tackle the approximation issue for the above-defined Donsker and Kac- 
Stroock processes by exhibiting a general convergence criterion (see Proposition I4.5|) , which 
will entail Theorem II .21 Eventually, we have added an appendix with material on fractional 
Sobolev spaces and the proof of a technical result needed in Section 13.41 

Remark 1.3. At first sight, the reader familiar with rough-paths type continuity results may 
be surprised at the absence of some 'Levy-area' term in our bound ([7]). Otherwise stated, the 
convergence of an approximation W"' towards W (with respect to some appropriate topology) 
is sufficient to guarantee the convergence of the associated solution. In fact, on this partic- 
ular point, the situation is very similar to the case of a one-dimensional SDE with so-called 
commuting vector fields, i.e., 

dYt = b{Yt) dt + (TiiXt) o dBl with (T-(x)ct,(2;) = a'Jx)ai{x) for ah i, j = 1, . . . , n. 

(10) 

It is a well-known fact (see for instance [55]) that under this commuting assumption, the 
solution Y of (jlOp appears as a continuous functional of the sole noise B (that is, no need for 
any Levy-area component). In a certain way. Equation ([2]) fits the above pattern. Indeed, for 
fixed X £ (0, 1), the noisy perturbation can be written as X^i^o[V^ei(a;)/(lf (x))] o d^l, i.e., 
we (morally) deal with n = oo and crj(-) = y/Xiei{x) f {•) in (fTOj) . So, at least at this heuristic 
level, our continuity result ([7]) becomes quite natural. In a more specific way, we will see that 
due to the commuting property, the Levy-area term arising from the rough-paths analysis of 
([2]) can be easily reduced to some continuous functional of W (Lemma 13. 8p . 

Remark 1.4. As we shall see it in Section [3l our proof of Theorem 11.11 heavilv relies on the 
properties of the fractional Sobolev spaces Ba,p, which we have recalled in the appendix. Un- 
fortunately, many of these properties become much more restrictive as soon as the underlying 
space dimension is larger than 2, as illustrated by the classical Sobolev embeddings. This 
accounts for our choice to stick to a one-dimension heat equation. Note however that our 
considerations on the theoretical study of (UD (Section [2D could be easily extended to a multi- 
dimensional setting. 

Remark 1.5. The results in this paper remain actually valid for any operator A of the form A = 
—dx{a-dx) + c, where c > and a : [0, 1] — )• M is a continuously differentiable function. Indeed, 
as explained in [10^ Section 2.1], such an operator A also generates an analytic semigroup of 
contractions and one can identify the domains T)[Ap) of its fractional powers with the spaces 
Ba^p, which is sufficient to follow the lines of our reasoning. 

Unless otherwise stated, any constant c or C appearing in our computations below is un- 
derstood as a generic constant which might change from line to line without further mention. 

2. The Stratonovich integral 
Recall that we are interested in the following mild equation: 

Yt = StiP + r St-u{f{Yu) o dWu), t G [0,r], iPeB, (11) 
Jo 

where {St)t>o denotes the strongly continuous semigroup of operators generated by —A with 
Dirichlet boundary conditions, and W is assumed to satisfy Hypothesis [TJ 

The integral appearing in (jlip is thus understood in some Stratonovich sense, an interpre- 
tation to be clarified in a convolutional setting, which is the main purpose of this first section. 
Once endowed with this interpretation, it turns out that ()lip reduces to a common mild Ito 
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equation with an additional drift term, and accordingly the existence and uniqueness of Y can 
be derived from well-known results (see Section r2.2|) . 

Note that the following regularity assumption on / will prevail throughout the section. 
Hypothesis 2. The function / : M — )• M zs bounded, of class and with bounded derivatives. 

2.1. The Stratonovich integral. In order to interpret J^St-uifiXu) ° dWu), we restrict 
our attention to a particular class of processes Y . Namely, we assume that, on some filtered 
probability space (O, J^, {J^t)t>o, P), {Yt, t £ [0, T]} is the unique valued mild solution of the 
following equation: 

dYt - AYtdt = V^dt + V^dW^ = e i3, (12) 

for some -F^-adapted random fields {V^*, t G [0, T]}, i = 1,2, with continuous paths in B (recall 
that B := L^{0, 1)). Moreover, we assume that 

supi?[||F,2|||] <+oo. (13) 

t<T 

In fact, such a process Y is explicitly given (see e.g. [9]) by 

Yt = St^l: + f St^sV}ds + /* St^s {V! ■ dWs) . (14) 
Jo Jo 

Now, a natural idea to define the integral in pi|) in some Stratonovich sense would be the 
following: introduce the kernel Gt-s{x,y) of St^s and, with the representation Q of in 
mind, set 

St^s ifiYs) o dW,) (x) " = " E ^ (/*^^*~^^^' *^^^^^^' ° "^^0 ' 

where the symbol * denotes the space variable and each integral j^{Gt-s{x, *)■ f{Ys), ej)is odpi 
is interpreted in the (classical) Stratonovich sense. Nevertheless, it is a well-known fact that 
the process Y defined by is not always a ;S-valued semimartingale (in other words, Y is 
not always a strong solution of (fH|) . see e.g. [9l Sec. 5.6]), making the definition of these 
integrals quite obscure at first sight. 

To overcome this difficulty, we consider a standard semimartingale approximation of Y: for 
every e > 0, let Y^ be the unique (strong) solution of 

dYl - A.Yfdt = Vldt + V^dWu Y^ = ^p, 

where 

-A, := i(Id-(Id-eA)-i) 

stands for the Yosida approximation of —A. In particular, — A^ defines a monotone and 
bounded operator which converges pointwise to —A (see e.g. [1]). Then, for every fixed e > 0, 
Y^ is a semimartingale, and we have (see e.g. O Proposition 7.5]) 

\\Tn suvE[\\Yl -YtWl] =0. (15) 
This extrinsic procedure will lead us to the following interpretation: 

Proposition 2.1. With the above notations, the family of Stratonovich integrals defined for 
allt£ [0,T],x G (0,1) by 

f St^s{f{Y!)odWs 

10 



(x):=liniVyA; / (Gt_,(x,*)./(y/),e,)eod/5n, (16) 



STRATONOVICH HEAT EQUATION 



7 



where the latter limit is considered in L'^{Q), converges in L^{Q]C^{[0,T];B)) as e tends to 0. 
Its limit, that we denote by S.-s {f{Ys) o dWg), satisfies the relation 

f St-s (/(K,) o dWs) = f St-s (/(n) • dWs) + f St-s (F/ . f'{Ys) ■ P) ds, (17) 
Jo Jo Jo 

where P{C) ■= \ Xlfcli ^k^kHY notation Sts (/(^s) • dWg) refers to the (usual) ltd 

integral. 

Thus, the Stratonovich integral in ()lip will henceforth be understood as in the latter propo- 
sition, in the class of processes Y satisfying an equation of the form ()12p . Note that the 
relation (|17p provides us with a familiar decomposition for the Stratonovich integral as the 
sum of an Ito integral and a trace term, and it must be compared with the decomposition for 
the (standard) Stratonovich integral. 

As a first step in the proof of Proposition 12.11 observe that the two terms in the right-hand 
side of PT]) are indeed well-defined processes in L'^{Q; B). This is a straightforward consequence 
of the boundedness of /, /', the trace-class assumption on W, and the fact that P defines a 
uniformly bounded function. 

We point out that, in the definition (jl6p . we first restrict the integral to {0,u) with u < t 
in order to avoid the singularity in the derivative of the kernel G. This will be clarified in the 
proof of the next lemma. 

Lemma 2.2. With the above notations, we have that, for all u G (0, t), 

JZV^j^J^ (G<„,(x,*) •/(y/),e,)i5 odf3ij (18) 

= / St-s{f{Y!)-dW,) {x)+ [St-s{nY!)-V^-P)]{x)ds. 
Jo J Jo 

Proof. For any fixed (u, x) E (0, t) x (0, 1) and j G N, the process s i— )• {Gt-s{x, *) • f{Yg), ej)e, 
s S [0, u], defines a (real- valued) semimartingale. Hence, we can use Ito's formula to assert 
that 



{Gt^six,*) ■ fiY!),ej)B = {Gtix,*) ■ f{^),ej)B + / {dtGtM^,*) ■ /(F/), e,-)^ 

Jo 

+ j\Gt-r{x,*){{^eYr' + V}) ■ f'{Y,') + ]^V? • /" (F/) } , e, )b dr 

OO 

+ y,V^k {Gt-r{x, *) • /'(y/) • V? ■ efe, e,)H d^^ (19) 
fc=i -^0 

The hypotheses on / and 1^*, and the fact that s < u < t, guarantee that all terms on the 
right-hand side above are well-defined. More precisely, using the spectral decomposition of G 
given by 

oo 

Gt-r{x,y) = Y.^~''"^''''~''^^k{x)eu{y), 

k=l 

one proves that, P-a.s., 

oo 



/ {dtGt-r{^,*) ■ f{Yn,ej)tsdr < C V (e" 
•^0 k=i ^ 
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and the latter is finite since s G (0, t). As far as the second pathwise integral on the right-hand 
side of ([IHD is concerned, we have, for instance, 



E 



f\Gt-r{x, *) ■ A,y/ • f'iY,'),e,)B dr 
Jo 



2 








< CE 


U'l 






Jo Jo 



Gt-r{x,ymeYr']{y)\^dydr 



<c(supi^[iiy/iii]) r^=L= 

^ r<T ^ Jo yt-r 



zdr < +00. 



Here, we have used the fact that < Gt~r{x,y) < (27r(t — r)) e ^(t-r) _ Similarly, one 
easily proves that the last term in ()19p is a well-defined square-integrable random variable. 

Plugging the expression ()19p in (jlSp and using the definition of the (standard) Stratonovich 
integral (see e.g. formula (3.9) in [191 P- 156]), we end up with 



M ,\Gt^s{x,*)-f{Y!),e,)Bodpl 




, ^ , (Gt_,(x,*)./(y/),e,)Bd/3i 





j=i J^ 



St-s {f{Y!) ■ dW, 
which concludes the proof. 

We can now go back to our main statement. 



(x)+ / [St.s{f'{Y^)-V!-P)\{x)ds, 
Jo 



□ 



Proof of Proposition Wl\ First, owing to the previous lemma, we have that the limit on the 
right-hand side of (fT6]) equals to 



St^s {f{Y!) ■ dWs) 



{x)+ / [St^s{f'{Y!)-V,'-P)]{x)ds. 
'0 



This can be proved using the bounded convergence theorem. Hence, the proof reduces to the 
two assertions: 



lim sup-E 

^->0 t<T 



f St-s (/(n') • dWs) - f St-s ifiYs. 
Jo Jo 



dWs 



0. 



(20) 



and 



lim sups' 



f St-sU'iY!) . ■ P)ds - f St-sU'iYs) ■ ■ P)d, 
Jo Jo 



0. (21) 
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Let US first deal with (|2U|) . By the isometry property of the stochastic integral, the boundedness 
of St^s and the assumptions on /, we have: 



E 



St^s if{Ys) ■ dWs) 



t oo 



Y,h\\St^s{[f{Y!) - f{Y,)]ek)\\lds 

k=l 

^ V r-T i-l 

<y^XkE / \f{Y'{s,y))-f{Y{s,y))\^\ek{vtdyds 
k=i I Jo Jo 

< Cj^ Afc supi^[||y/ - YtWl] < C supE[\\Yf - Yt\\l] , 



k=l 

^00 



t<T 



t<T 



upon recalling that X^^i < oo. 

In order to prove ()2ip . we use the Sobolev embedding L^(0, 1) C B_i_^ and the assumptions 

4 

on / and V"^. In fact, we have 
ft 

E 



[ St-s{[f'{Y,^)-f'{Ys)]-V,'-P)ds 
Jo 

< C\\PU^ f\t-sn-^E[\\[f{Yl)-f'{Ys)]-V,^LAds 

Jo 

< C\\PU^ f \t- sn-^ EiWifiYl) - f'{YMB\\V,'\\B]ds 

Jo 



(22) 



< c (sup£;[||y/ - y,|||]V2) (supEmYB]' 



t<T 



1/2 



t<T 



(recall that || • refers to the supremum norm on [0, 1]). Therefore, by the assumptions 
on V'^, the convergence (fT5|) guarantees that ([20]) and (f2T|) hold, and this lets us conclude the 
proof. □ 

2.2. Existence and uniqueness of solution. With the notations of the previous section, 
consider the following mild (Ito) equation: 

Yt = StiP + f St^s {f{Ys) ■ dWs) + f St-s{f{Ys) ■ f{Ys) ■ P) ds, (23) 
Jo Jo 

where we recall that -P(^) = YlT=i ^k^kii)"^ ■ Hypotheses [1] and [2] allow us to apply standard 
methods and guarantee that this equation admits a unique L^(0; B)-valued solution Y (see [U]). 
In particular, we observe that Y solves an equation of the form (|12p with = f'iYt) ■ fiXt) ■ P 
and = fiYt) and these random fields fulfill the assumptions specified in the previous section. 
Thus, for all t G [0, T], we can define the Stratonovich integral Jq St-s (fiYg) o dWg) through 
Proposition 12.11 and we know that 

f St^s{f{Ys)odWs)= f St^s{f{Ys)-dWs)+ f St^sif'iYs)- f{Ys)-P)ds, 
Jo Jo Jo 

which yields that Y is also a solution of (|lip . 

Conversely, due to p!7|) , it is readily checked that any solution of ([TT]) in the class of processes 
satisfying an equation of the form ()12p is also a solution of ()23p (use the uniqueness of V^, V'^ 
in ()14p ). This provides us with the following existence and uniqueness result. 
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Theorem 2.3. Assume that Hypotheses [J\ and\^ are both satisfied and that tp B. Then, 
there exists a unique B-valued process \Yt, t G [0,T]} which solves 

Yt = Stij + f St^u{f{Yu) o dWu), t G [o,r]. 
Jo 

Moreover, Y has a version with continuous paths and it holds that supj<y < oo. 

3. A ROUGH-PATHS TYPE ANALYSIS OF THE EQUATION 

Let us now turn to the proof of Theorem ll.il As announced m the Introduction, our strategy 
is based on a rough-paths type expansion of the equation. Accordingly, a few ingredients taken 
from the so-cahed convolutional rough paths theory, that is rough paths theory adapted to mild 
evolution equation, must be introduced in the first place. 

3.1. Tools from (convolutional) rough paths theory. We gather here some preliminary 
material borrowed from [16] (see also |1H llOj). As underlined in the latter references, a key 
point towards a fruitful pathwise analysis of dH lies in the following elementary observation: 
due to the semigroup property St+t> = St ■ Sf, it holds that, for any s < t, 

Yt-Ys = f St-uU{Yu) o dWu) + r [St-u - Ss-u] U{Yu) o dWu) 

J s Jo 

t 

St-u ifiYu) ° dWu) + atsYs, where ats := Sts - Id . 

Otherwise stated, by setting {5Y)ts ■= (Yt — Yg) — atsYg, the equation (jH can be equivalently 
written in the convenient form: 

yo = V', {SY)ts= f St-u ifiYu) odWu), 0<s<t<l. (24) 

J s 

This should be compared with the behaviour of solutions to standard (stochastic) differential 
equations: if Xt = a + aiXu) dBu, then i6X)ts '■= Xt — Xg = aiXu) dBu- Then, in 
a rough-paths setting, we are naturally led to extend the definition of 5 to processes with 2 
variables, as follows: 

Notation 3.1. For all processes y : [0,T] — )• ;B and z : B, where S2 '■= {(s,i) £ [0,r]^ : 

s < t} denotes the two-dimensional simplex, we set, for s < u < t & [0, 7"]- 

iSy)ts ■■= yt - ys, iSy)ts ■= iSy)ts - atsys = yt- Stsys, (25) 

i^z)tus '■= Zts — Ztu — St-uZus- (26) 

To make the notations (|25 p - (|26p even more legitimate in this convolutional context, let us 
point out the following algebraic properties: 

Proposition 3.2. For any y : [0,T] — )• B, it holds: 

(i) Telescopic sum: 5i5y)tus = and i5y)ts = YhZq St-t,+Ah)u+iU for any partition 
{s = to < ti < . . . < tn = t} of an interval [s, t] of [0, T] . 

(ii) Chasles relation: if Jts '■= jl St-u [yu ■ dWu), then 6 J = 0. 

Both points (i) and (ii) are straighforward consequences of the semigroup property. Now, in 
accordance with the new expression (j24p for the equation, a (5-version of the classical Holder 
norm must come into the picture. To this end, fix a subinterval / C [0, T] and a Banach space 
V. Then, if y : / ^ F and A > 0, set 

N[y;e\T,V)]:= sup (27) 

s<t&I \t — s\ 
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and define C'^(/; V) as the set of processes y : I ^ V such that M[y; C^{I; V)] < oo. 

As we will see it in the sequel, a proper control for the expansion of St-u{f iXu) ° dWy) 
also requires the extension of both definitions ^ and (j27p to processes with 2 or 3 variables. 
Precisely, \i z : S2 and /i : ^3 — )• F, where ^3 := {(t, u, s) G [0, T]^ : s < it < i}, we set 

AA[z;C2"(/;y)]:= sup N[h-CUl;V)]:= sup (28) 

s<tei \t — s\ s<u<tei \t — s\ 

and we define C2 (/; V) (resp. C^{I', V)) along the same lines as C^{I; V). Observe for instance 
that if ?/ G €2(1; C{V,W)) and z G C2{I;V), then the process h defined as htus = UtuZus 
{s<u<tel) belongs to C^^'^{I; W). 

Note that when / = [0, T], we will more simply write Cj^{V) := Cj^{I;V) for k £ {1,2,3}. 
Besides, from now on, we use the following convenient notation for products of processes. 

Notation 3.3. If g : Sn ^ C{V,W) and h : Sm — ^ W (with n,m £ {l,2,3}j, we define the 
product gh : Sn+m-i W by the formula 

With this convention, it is readily checked that g : S2 ^ C{B^,Ba) and h : [0,T] — )• iS^, 
then 6{gh) : S3 ^ Ba obeys the rule: 

5igh) = i5g)h-g{5h). (30) 

To end up with this toolbox, let us report what may be seen as the cornerstone result of the 
convolutional rough paths theory, namely the existence of (some kind of) an inverse operator 
for 6, denoted by A, and which will play a prominent role in our forthcoming decomposition 
(pI2|) . In brief, this operator allows us to get both a nice expression and a sharp estimate for 
the regular terms, i.e., the terms with Holder regularity strictly larger than 1, that arise from 
the expansion of St-u{f{yu) ° dWu) (see in particular the proof of Lemma [3. lip . 

Theorem 3.4. Fix an interval I C [0,T], a parameter k > and let /i > 1. For any 

h £ C3 (/; B) n /to (5, there exists a unique element 

A/iGn,g[o,^)Cr"(I;e,) 

such that 6{Ah) = h. Moreover, Ah satisfies the following contraction property: for all a G 
[0,/i), 

AA[A/i;Cr"(/;Sa)] < Ca,^,M[h;C^{I;B)]. (31) 



The proof of this result can be found in [16\ Theorem 3.5]. 

3.2. A rough-paths type expansion of the solution. We are now ready to settle our 
reasoning, which applies to a smooth enough vector field /: 

Hypothesis 3. The function / : M — )• M in |^]) is of class , hounded and with hounded 
derivatives. 

Our main task will actually consist in establishing the following pathwise decomposition for 
the solution Y to (IMD: 

Theorem 3.5. Assume that hoth Hypotheses\^and\^hold true. Fix ^ G (^, ^ + r]) and assume 
that ip £ B-y. Then the 5 -variations of the solution Y to \24\ ) can he expanded as 

{5Y)t, = f St^u{f{Yu) o dWu) = LlUiYs)) + Lf^'^ifiYs) ■ fiY,)) + Ats{R^), (32) 
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where we have set, for all s < u < t, 



lY{^):= I St-u{^-dWu), (33) 



^ (^) := / St-u ■ {SW)us ■ dWu) + / St-u ■ P) du, (34) 

J s J s 

and 

' tus 



Rl, := -6{L^f{Y) + L^^(/(y) • /'(F))),,,. (35) 



The theorem must be read as follows: in the expansion of St~u{fiXu) ° dWu), we can 
exhibit a main term, namely 

Lir(/(ys))+Lir^(/(n)-m)), 

and a residual term Ats{R^) with Holder regularity strictly larger than 1, in the sense of 
Theorem 13.41 (take a = in (j3ip ). Besides, from the decomposition (j32p . we can somehow 
conclude that the whole dynamics induced by W is "encoded" through the two (stochastic) 
operator- valued processes and . So, before we turn to the proof of ()32p . let us 

elaborate on the properties of these two processes. 

3.3. The couple {L^ , L^^). At this point, we consider and as stochastic linear 

operators acting on the space of smooth functions ip. The following (straightforward) relation 
accounts for the algebraic behaviour of the couple {L^ , L^^): it is the convolutional analog 
of the classical Chen's relation between a process and its Levy area (see |15|). 

Proposition 3.6. The processes and jj^^ obey the following algebraic rules: For all 
s < u < t and all smooth function (p, 

{6L'^)tus{^)=0 , {6L'^'^)tus{^) = LZ{y^-{6W)u,). (36) 

Now, it matters to identify the regularity properties of and as 2- variables processes. 

A first clue in this direction is given by the following (a.s.) regularity result for the noise W 
itself. 

Lemma 3.7. Under HypothesisUl one has (a.s.) W £ C^~^{Bri^2p) for every integer p > 1 
and every small e > 0. 

Proof. By using our forthcoming Proposition 14.31 we deduce that 

^[iiWi.iiSy <c^p,.^[i('5/?)*.i'"l <t^p,.ii-^r, 

for all q > 1, and the result is now a straightforward consequence of the Garsia-Rodemich- 
Rumsey Lemma l6. II (take 5* = 6 and R = 6W in the latter statement). □ 

Our second ingredient towards the regularity properties of {L^ , L^^) relies on two suc- 
cessive observations. First, due to their relative simplicity, the two expressions ([33|) -([34 |) can 
be integrated by parts. Then, owing to the some obvious commuting properties, we can turn 
j^ww g^j^ easy-to-handle functional of 6W. This is what we propose to detail in the proof 
of the following Lemma. 

Lemma 3.8. For every smooth function if and all s < t £ [0,T], the following formulas hold 
true (a.s.): 

L^i^) = St-siip • {6W)ts) - ^ ASt-u{^ • {6W)tu) du, (37) 
L'^'^iv) = \ [St^si^ ■ {5W)l) - J' ASt^u ■ [{SW)l + 2{6W)tu ■ {SWU]) du^ . (38) 
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Proof. With the expansion ([3]) in mind, it is easily checked, by setting := X^^i ^/XiPlci 
that Ljf ((^) = hniAT^^oo Ei=i L ^t-uiv ■ dW^) and 
ft ft 



St^ui^ ■ {SW)us ■ dWu) = hm / St-ui^-idW^'U-dW^), 

where the Hmits are taken in L^(r2, B). The proof then reduces to apphcations of Ito's formula 
and we only elaborate on ()38p . For fixed i, j G {1, . . . , A^}, apply Ito's formula to the (random) 
function Flj : [s,t] xM. xM. ^ B defined by 

so as to deduce 

+ / St-u{<f-ei-ej){5f3'')usdl3i+ / 5't„„(v3- • ej)((5/3-^')„s + / St-u{f-ei-ej)du. 

Js Js Js 

By taking the sum over we deduce the formula 
St^siv ■ - ASt^u{^ ■ [(<5W^)L - (<5Vr^)L]) du 

ft ft N 

= 2 St-u{ip-i6W'')us-dW^)+ / St-u{^-{^Xie])) 

J S J S .-I 



i=l 



and by passing to the limit (in L {^},B)), we get 



^JT'^M = \[st-s{^ ■ {6W)l) + I ASt-u{y^ ■ [{SW)l - {5W)l]) dn\. 
Formula ([55|) immediately follows. 



□ 



We are now in a position to extend both and to larger classes of functions ip and 

retrieve the following (a.s.) bounds, which will be at the core of our identification procedure: 

Proposition 3.9. Under the hypotheses of Theorem \3.5\ for any small e > 0, there exists 
e > and p > 1 such that (almost surely) 

M[L^;Cr\C{B, B))] + AA[L^; 4^-^+''^"'(/:(fii/2. ^7))] + Ml^^ ■,Cr\C{B^, B^))] 

<Ce,e,pM[W-d~'{B^,2p)l (39) 



M[L'^'^;Cl^'{C{B,B))] +AA[L^'^;C('-^+'')-^(/:(^i/2,fi-,))] +AA[L^^;C| \C{B^,B^))] 

<Ce,e,pM[W;d-'{Br,,2p)]\ (40) 

for some constant Ce^e,p- 

Note here how important the assumption 7 < ^ + 77 in Theorem 13.51 to ensure that (^ — 7 + 
r/) — e > for any small enough e > 0. 



Proof. In fact, thanks to the representation formulas (|37p - ()38p and the pathwise regularity 
of W f Lemma 13. 7p . all of these bounds can be derived from the classical properties of the 
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fractional Sobolev spaces (see Appendix A). For instance, owing to (|72|) . one has, for any p > 1 
and a > ^, 

\\^ • {SW)ts\\l3-c < Ca,p\\(.^W)ts\\L^P{o,l)M\B, (41) 
so that for any i small enough, 

\\St-s{^ ■ {6W)ts)\\B < Ca,p\t-s\-''\\{6W)ts\\L^P(o,i)\MB (use ([TO])) 
In the same way, 

WASt^uiif ■ iSW)tu\\B < c^,p\t-s\-^-''\\{SW)ts\\ l2p(o,i)IIv'IIb (use (HDD) 

< cc.,p,e\t - s\^i~'"'''>~^ M[W;C-2-%Br,,2p)] M\b- 
By taking a small enough, i.e., p large enough, we get the expected bound, namely 

M[L'^;Cr'{C{B,B))] < Ce,e,pM[W-d-'iB.,,2p)]- 

The other estimates for can be proved along the same lines. As far as is concerned, 

observe for instance that if e > is small enough, then one has 

\\ASt-u{^ ■ {6W)tu ■ {dW)us)\\B-, 

< c\t- nl^^''-^)-^ y ■ {6W)tu ■ {SWUWb, 

< c\t - 4"^-^^-^ \\{6W)tu ■ i6W)us\\BM\B,^, {nsem) 

< c|t-n|(''-^)"i||(5t^)j„||e,,J|(5T^)ns||i3,,4ll<^llBi/2 (use (Ell)) 

< c. It - u\^-2+^-^-'^-' \u - s\"^-'M[W;d-%B,,^)fy\\B,,„ 
which entails that 

AA[L^^;4^-^+'')-^(£(Si/2,e^))] < Ce,e,p^f[W■d-%B,,,)]''. 
The (analogous) proofs for the other bounds are left to the reader. □ 

3.4. Proof of Theorem 13.51 First, we need to justify that the right-hand side of the decom- 
position (|32|) is well-defined. This will rely (among others) on the following a priori controls 
for the solution Y. For the sake of clarity, we have postponed the proof of this statement to 
Appendix B. 

Lemma 3.10. Under the hypotheses of Theorem \3. 51 one has (almost surely) 

Y eC^'^{B^)nC\B-,), (42) 

:= 5Y - L^f{Y) G C|^''(S). (43) 

Recall that according to our convention (|29p , the definition of in (|43|) must be understood 
as := {SY)ts - L'^ifiYs)) for every s<te [0,T]. 

Lemma 3.11. Under the hypotheses of Theorem \3.5l let Z be the process given hy Zq = and 

{5Z\s = LYAfiYs)) + ^ir^(/(n) • f'iY,)) + ku{B^). (44) 

Then, almost surely, Z is well-defined as an element of C"^^ (Boo) riC^{B-y), and there exists a 
constant A > such that for any subinterval I = [^1,^2] C [0,r], one has 

M[Z;C^^iT,Bo.)] +M[Z;C%T,B^)] < \\Z,,\\b, + cwM^'MiY; Qil)], (45) 

where we have set 

N[Y- Q{I)\ := Af[Y; C'^{T, B^)] + AA[F; C°(I; B^)] + AA[K^; C|+'(/; B)] 
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Proof. First, according to Theorem \'6A\ we need to justify that € Ct^i^) ^ov some > 1. 
To this end, expand R using the algebraic rules (j3U|) and (|42|) . which gives 

RL = lZNus + LZ'^SifiY) • f'{Y)U, (46) 

with Nus ■■= S{f{Y))us - {6W)us ■ f{Ys) • f'{Ys). Thanks to (gOD and (gS]), it is readily checked 
that L^^6{f {Y) ■ f'{Y)) G Cl'^'^'^~%B) for any small e > 0, since 

\m{Y)- f'{Y))us\\B < c\\{6Y)us\\b < c{\\{5Y)us\\B + \\ausYs\\B} 

< c{\u - sf''^Af[Y;C^\Boo)] + \u- s| W[y; C^S^)]}, 

where we have used ([7T]) to get the last inequality (recall that ats '■= Sts — Id). 

Then, as far as N is concerned, let us expand N using standard differential calculus, which 
provides us with the expression 



dr f'{Ys + t{5YU) ■ {Kl, + ausYs + K^ifiY,))} 

Jo 

+ f dr [f'{Ys + r{5Y)us) - /'(n)] • m)us ■ f{Ys), (47) 
Jo 

where the additional operator-valued process L""^ is defined by 

Lt^i^) ■■= f atu{^ ■ dWu) = ats{^ ■ {6W)ts) - f ASt-u{v ■ {SW)tu)du. 

J s J s 

-—e -+e 

Now, since G q {C{B,B)), it is sufficient to prove that N £ Cf iB) for some small 
e > 0. But, with the expansion (j47p in hand, this becomes an easy consequence of the a priori 
controls given by Lemma 13.101 together with the regularity property: 

AA[L'^^;C^^+''^~'(£(fii/2,e))] < c,,,,,J^[W;d-'{B,,2p)], 
derived from ()73p . Note in particular how important the assumption 7 > 1/2, insofar as, by 

m, 

\\f'iYs + ri6Y)us) ■ {ausYs)\\B < CjIKsY^Wb < C\u - AA[y; C°(^^)]. 

We are thus in a position to apply A to , and so Z is properly defined through (|44p . The 
regularity of Z and the bound (|45p are immediate consequences of (j39p - (j40p and the contraction 
property ([3T]) of A. The details are left to the reader. 

□ 

Remark 3.12. Although not optimal, the two regularity results (j42p and (|43p are thus sufficient 
for us to prove that the right-hand side of the decomposition (j32p is indeed well-defined. We 
also retrieve an important stability phenomenon here: Y and Z both belong to the same space 
C'^^{Boo)riC^{B.y). A posteriori, this accounts for our choice in favor of this particular topology. 



We can eventually proceed to prove Theorem 13.51 

Proof of Theorem 13.51 We need to identify the increments of Y with those of the process Z 
defined in Lemma 13.111 To do so, we naturally rely on some expansion of the right-hand side 
of ([23|) . Precisely, we have that 

r St-u{f{Yu) ■ dWu) + f St-u{P ■ f{Yu) ■ f'{Yu)) du 

J s J s 

= L^ifiY,)) + Ljr^(/(y.) • /'(n)) + jYs, 
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with 

jZ:= [ St^u{P-S{f{Y).f'{Y))us)du+[ St^u{Nl ■ dWu) , (48) 

J s J s 

where the process = d{f{Y))ts — {6W)ts ■ f{Ys) ■ f'{Ys) has already been considered in the 
proof of Lemma 13.111 Therefore, with this notation, it holds that 

6{Z-Y) = AUR'')-jrs. 

Now, by the contraction property (jSip . we know that A{R^) G Cl^^ (B) for some fii > 1. Besides, 
with the same ingredients as in the proof of Lemma 13.101 (Burkholder-Davis-Gundy inequality 
plus Lemma 16. H see Appendix B), we can easily lean on the expansion (j47p of N to prove 
that G Cf (S) for some fi2 > 1 (note that 5J^ = R^). Consequently, 5{Z -Y) e C^{B) 
with fi = inf(^i,/i2) > and this entails that 6{Z — Y) = 0. Indeed, for any partition 
= {s = ti < . . . < tn = t} [s,t], one has, due to the telescopic sum property reported 
in Proposition \3.2\ 

\\d{z - Y)ts\\B < \\^(^ - Y)u+iuh < - ur < i^mT"' I* - ' 

i i 

and we conclude by letting the mesh |'P[s,t]| := maxj — ti\ tend to 0. □ 

As a straightforward consequence of the decomposition ()32p , we can exhibit an almost sure 
bound for Y in terms of W . Indeed, by plugging the estimate (|45p back into the equation, we 
deduce that for any subinterval / = [^1,^2] C [0,r], 

M[Y-C^\I-Boo)] +N[Y-C\I-B.,)] < \\Y,,U,+Cw\I\''M[Y-Q{I)] 

for some constant A > 0, and similar estimates for = [f {Y) ■ f'{Y)) + A{R^) finally 
show that 

AA[y;Q(/)] < \\Ye,\\B,+Cw\lMY;Q{I)]. 
At this point, a basic patching argument easily leads us to the following statement: 

Corollary 3.13. Under the hypotheses of Theorem lS.Sl there exist e > and p > 1 such that 

M[Y;Q{[0,T])] < G,,p{\ms„M[W;d~%B,,2p)]) (49) 

for some deterministic function Gg^p : (M"*")^ — >■ bounded on bounded sets. 



3.5. Comparison with smooth solutions. The previous considerations will allow us to 
prove our continuity result (Theorem II. ip and for this purpose, we first go back to the case 
where the driving noise is an absolutely continuous process W (with values in ^r;,2p); assumingly 
defined on the same probability space as W . In this situation, our mild equation is naturally 
understood in a pathwise sense as a classical (Riemann-Lebesgue) mild equation, i.e., 

Yt = Sti^ + f St-u{f{Yu) ■ dWu) = St^ + /* St-u{f{Yu) ■ Wi) du, (50) 
Jo Jo 

and the (pathwise) existence and uniqueness of the solution Y follows from standard PDE 
results. The key step towards a comparison between Y and Y lies in the following result, 
which points out the similarity between the couple {L^,L^^) at the core of the previous 
considerations and the couple (L^ , L^^) constructed from W: 
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Lemma 3.14. Define the operator-valued processes and in the classical Riemann- 

Lehesgue sense as 

lS{^):= 1^ St-u{^-dWu) , lS'^{^):= I' St-u[^-{6W)us-dWu), (51) 

for every smooth function ip. Then both formulas ( |^7| j and ll'^8\) remain valid when substituting 
W for W , and accordingly the bounds 139\) and hold true for W as well. 

Proof. It suffices to replace the use of Ito's formula in the proof of Lemma [3.81 with standard in- 
tegration by parts. Indeed, as an absolutely continuous process, W obeys the rules of standard 
differential calculus and one has for instance 

:= l\sW)us ■ dWu = \{6W)l. 

Consequently, it holds that 



i: 



St-ui^ ■ i5W)i) - - / ASt^,,iip ■ [i6W)l - i6W)l])du 



= ^St-u{y^ • {dW)l) - ^ 

= ^St.^{^ ■ {6W)i) - ^ £ ASt-u{^ ■ [{6W)l + 2{6W)tu ■ {SWU] ) dn, 

which precisely fits the pattern of (|38|) . □ 

Another consequence of the similarity between {L^ , L^^) and (L^, L^^) through the two 
formulas (|37p and (|38p is a set of (readily-checked) Lipschitz-type bounds: with the notations 
of Proposition 13. 9| one has, for some polynomial expression c^^y, 

MIL"^ - L^;Cr\C{B,B))] < c^^M[W - Vr;C^^'(^,,2p)], (52) 

AA[L^^ - L^^;Ci-^(£(^,e))] < c^^N[W - W;d-%B,,2p)], (53) 
and this bound remains valid for all of the other topologies involved in Proposition 13.91 
Then, as far as the solution Y is concerned, note that 

(Sns = Ltf{Y,) + Ljf ^(/(y,) • f'{Y,)) + jI 

with := tsSt~n{[5f{Y)us-{5W)us-f{Ysyf'{Ys)]-dWu), audit is obvious in this (absolutely 
continuous) situation that S (■S) for some /i > 1. Therefore, we can easily follow the 



lines of our previous identification procedure (see the proofs of Lemma 13.111 and Theorem 13.5 
in order to exhibit a similar formula for the J-variations of Y: 

Lemma 3.15. Under the hypotheses of Theorem assume that ip G B^. Then the 6- 
variations of the solution Y to Ii50\) can be expanded as 

{6Y)t, = Lt{f{Ys)) + Lt^{f{%) ■ f'iYs)) + At.(#), (54) 

where i?L := -6{L^fiY) +L^^{f{Y) ■ f'{Y)))^^^. In particular, the bound remains 
valid for Y when replacing ip (resp. W) with ip (resp. W). 
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With these identifications in hand, the proof of Theorem 11.11 becomes a matter of a standard 
rough-paths argument, and we only sketch out the main steps of the procedure (see e.g. the 
proof of Lemma 5.2] for further details on the computations). 

Proof of Theorem \l.l\ In order to compare Y with we can now rely on their respective 
decompositions ([32]) and ([51]) . By setting g := ff, we get that 

6iY - Y),, = { [Lir - Lt]f{Y,) + [LJT^ - if^] {g{Y,)] 

+ {LJT [/(y.) - fiYs)] + lS'^ [a{Y,) - g{%)\ } + hsiR" - (55) 

with a similar splitting for — BX (based on the expansion ()47|) ). Now, as in Lemma 13.111 
we consider the following appropriate topology: 

N[Y - y; Q(/)] := M[Y - Y; C^\I; B^)] + M[Y - Y; C\l; B^)] + J\f[K^ - K^;Cl''\l; B)]. 

By using the decomposition ([55]) and the bounds ([52]) - ([53]) . standard differential calculus shows 
that for any subinterval / = [^1,^2] of [0,T], 

M[Y-Y-QiI)] 

< Cw,w,^,^ {ll^^i - ^iih, +^[W - W-d"%B,,2p)] + \I\'M[Y - Y; Q(/)]} , 

for some constant A > 0. As in Corollary I3.13| we can then rely on an elementary patching 
argument to reach the global bound ([7]). 

□ 

Remark 3.16. The above strategy sheds new light on the classical Ito-Stratonovich correc- 
tion phenomenon arising in the approximation of stochastic heat equations. Indeed, on the 
one hand, it emphasizes that the convergence of Y towards Y reduces to the convergence of 
(L^,L^^) towards {L^ , L^^), and on the other, continuous bounds such as ([53]) clearly 
highlight the relevance of the Stratonovich interpretation of in this context. In a way, 

the correction phenomenon is therefore more directly observed through the decomposition (|34p 
of as the sum of an Ito integral and a trace term. 

4. Approximations in law 

We now aim to prove our approximation result, that is Theorem 11.21 Thus, from now on, 
we assume that the hypotheses in Theorem 1 1 . 2 1 are all satisfied. Recall that the approximation 
processes involved in this statement, namely the Donsker and the Kac-Stroock approximations, 
have been specified in the Introduction (see ([8]) and dS])), as well as the notations W and 
Besides, in this part of the paper we take T = 1 for the sake of simplicity. 

4.1. Preliminary results. As a first step towards Theorem ll.2| we need to check that the 
processes we have constructed via W are indeed well-defined. To do so, we will make use of 
the following bound. 

(n) in) 

Lemma 4.1. Fix n > 1. Let X\ ,...,Xn be independent centered random variables with 
moments of any order and fn : {1, . . . ,n} — t- M. Then for every r > 1, there exists a constant 
Cr which only depends on r such that 



E 



i-i J V ,_i / V l<«<n 



('^)|2rj 
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This inequality can be easily deduced from the following result, which is clear for r = 1 and 
was proved by Rosenthal for r > 1 (see |25t Thm. 3]). 

Theorem 4.2. Let Yi, . . . ,Yn be independent centered random variables satisfying < 
oo, where r > 1. Then, there exists a constant Cr such that 



E 



n -1 n / " \ 

<Crma^{Y,E\Y,r,(Y,E\YA }■ 

i=l ^ i=l ^ i=l ^ 



The transition from real-valued to S^.2p-valued processes will be ensured by the following 
result. 

Proposition 4.3. Let (Xfc)fc>i he a sequence of centered i.i.d. random variables on some 
probability space {Q,J^,P). Assume that each Xi has moments of any order, and consider a 
sequence {Xk)k>i of positive numbers such that 'Ylik>i < oo for some (fixed) r] > 0. Then, 

for every p,q > 1, the random series of functions \/AfcX,fcefc converges in L'^P'^{^1, 13ri^2p) to 
an element X which satisfies 

E[mX'J<Cp,g,x,vE[\X,\'P'^], (56) 
for some constant Cp,q^\^^ which only depends on p, q and ^fc>x k^^ . 

Proof Set X" := J2k=i Vh^^kek and observe that ||X™-X"||b^ = ||X('"''^)'''||i2p(o,i), where 
we have set x('^'")'''(0 := YJk=n+i k^"^ V^^XkCkiO- Then, by Jensen's inequality, 



1 \ <?■ 



< /^d^^[|x(™'")'^(e)P^"'] 
Jo 



and thanks to Lemma l4.H we get 

„i . m 

^[||X(-,n,),,||2p. 1 < C,,,E[\X,\'r'^] dd E ^kk'^eUC? 

^k=n+l 

< Cp,,E[\X,fP'i]( ^kk'^] (57) 

^ k=n+l ^ 

due to the uniform bound ||efc||e^ < \f2. In particular, i?[||X'" — X^H^^*^^^] tends to zero as 

both m and n tend to infinity, so that X" converges in L'^P'^(il,Br]^2p)- The bound (|55|) can of 
course be derived from (|57|) . □ 

In particular, due to Hypothesis [H we can conclude that = W(S'"''') and W"' = W(0"' ) 
are indeed well-defined processes with values in Bri^2p- Let us now get a little bit closer to 
the assumptions of Theorem 11.11 by checking that in both cases, admits an absolutely 
continuous version. 

Lemma 4.4. For any fixed n > 1, both the Donsker approximation = W((S'"''') and the 
Kac-Stroock approximation = W(0"' ) have an absolutely continuous version with values 
in Brj^2p, for all p > 1. 

Proof. Since the (deterministic) approximation grid for S^'^ does not depend on k, it is easily 
seen that 



W(5"'-)t = W(5'^'-)± (i- -) • {W(S"'-)i+i - W(5"'-)±} ifiG + l 



-n n 

In particular, W(5'"' ) is a piecewise linear process (with values in Bri^2p) and accordingly it is 
absolutely continuous. 
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As far as the Kac-Stroock approximation is concerned, first we can see that it has a contin- 
uous version with values in i3^,2p- Indeed, applying Proposition] 



E 



|<5(W(r'-))ts 



< CE 
= CE 



t 



,2pq 



2pqi 



n(-l)^+^(™)(in <CnP'^\t-s 



2pq 



For the sake of clarity, we will also denote by W(^"'') this continuous version. To prove the 
existence of an absolutely continuous version, we will see that with probability 1, 

W(r' )t = / W(^"'% ds, for any t G [0, 1], (58) 
^0 

where 9^ := ^/n ■ (— l)C+^("*). Indeed, thanks to Proposition 14. 3t W(^"'')t is well-defined for 
every t G [0, 1] as an element of L^^(r2, 5^^2p) and 



E 



|W( 



< 



E[\\-W{9"'-) 



n,-\ u2p 
sWb. 



V.2p- 



2p 



ds < Cr, Elm^p]]'" ds < 



oo. 



As a consequence W(0"''') is (a.s.) Bochner-integrable. Moreover, for each t G [0, 1], 



N ft N 

W(r' )t= lim Vv^^,"''=efc= lim / ( V yA^C'efc)^^ in ^^^(O, ^^,2p) 
k=\ ■'^ k=l 



and 
E 



N 



n,k 

k 



'ek)ds- I W{9"'-)sds 



k=l 







Br],2p 



^ 7. — AT I 1 



|2p 



k=N+l 



Jri, 2p 



ds 



<CJ sup E[\9^\'P'>])( ^kk^"" 



sg[0,l] 



k=N+l 



by similar arguments as in the proof of Proposition 14. 3i Since the last expression tends to 
as — 7- oo, we obtain that for each t G [0, 1] 



W( 



ds a.s. 







Thus, since {W(6'"'')j, t G [0, 1]} and {/g W(^"'% ds, t G [0, 1]} are both continuous processes, 
we can conclude that 

ft 



P{w{9^'-)t = ^ 



W(9"' )sds, yt G 



[0,1]} 



1. 



□ 



4.2. A general convergence criterion. One of our key ingredients to prove Theorem 11.21 
via Theorem 1 1.1 1 lies in the following statement, which puts forward sufficient conditions for an 
approximation of the noise (defined on the same probability space) to converge with respect 
to the topology involved in d?]). 

Proposition 4.5. Let (/3")„>i be a sequence of centered processes and /3 a Brownian motion, 
all defined on a same probability space {^l,J^,P), and such that the following two conditions 
are satisfied: 

(i) For every integer p > 1, there exists a constant Cp such that for all s,t £ [0, 1] and all 
n>l, 

E[\P^-P^fP]<Cp\t-s\r 
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(ii) For every integer p > 1, there exists a constant Cp such that for all n > 1, 

sup E[\I3^ - Ptl^P] <Cpn-''P, 

te[o,i] 

for some fixed parameter u > 0. 

Then if we consider independent copies {I3'^'^)k>i (resp. (/?^)/c>ij o/ /3" (resp. (3) on a same 
probability space, we have that, for any integer p > 1 and any e > 0, 



W(/3 );C2--(^,,2p)] 

n— >-oo 



a.s. 



Let us first see how to combine the above conditions (i) and (ii) so as to exhibit convergent 
bounds in Holder topology. 



Lemma 4.6. Under the hypotheses of Proposition \4-5l for all integers n,p > 1, all e G (0, 1) 

and s < t & [0, 1], one has 

for some constant Cp which only depends on p. 

Proof, li \ t — s\ < n~'^, then due to the condition (i), it holds that 

E[\6{f3'' - f3)tsfP] < Cp{E[\^^ - f3^\^P] +E[f]t-f3s\^^]}<Cp\t-s\P<C, 



t 



On the other hand, ii \t — s\ > n one has, thanks to the condition (ii). 



E[\6{f3'' - (3U^P] < Cp sup E[\I3^ - < CpU-'^P < Cp- 

te[o,i] 



□ 



Proof of Proposition \4-5\ By using successively Proposition 14.31 and Lemma 14. 6t we get, for 
any q>l, 

1 r 1 \t — q\(^~^)pi 

|5(W(/3"'-) - W(/3-))^J||^'^ p lAr«" _ |2pg — *i 



E 



^ Cp^q^n E 



|5(/3"-/3)t. 



< a 



^pqve 



We are thus in a position to apply the Garsia-Rodemich-Rumsey Lemma 16.11 (with 6* = 5) 
and assert that, for q large enough. 



E 



■aa[w(/3"'-) - winc'^-'iK^p)?"' 



< c 



p,q,v 



— Cp^q n IT' 



E 



||<5(W(/5"0-W(/3-)),J||7,^ 



[0,T]2 
epqu 



\t - s|2P9(|-£)+2 

\t - sY"^'-"^ dsdt < C 



dsdt 



[0,T]2 



n 



-spqv 



As a result, it holds that 

P(AA[W(/3-'-) - W(/3-);C^^(^,,2p)] > n— /^) < C^,,,, n'^^'^'^/^, 
which, thanks to the Borell-Cantelli Lemma, leads us to the conclusion, that is 

AA[W(/3"'-) - W(/3-);C^-^(i3^,2p)] ^ a.s. 

as n tends to infinity. 



□ 
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Example: As an immediate illustration of Proposition 14. 5t let us consider here the Wong- 
Zakai approximation of a given noise W satisfying Hypothesis [TJ Precisely, set 

:= W± + n ■ (t - -) ■ {W^+l - W^} forte 
and denote by the solution of the equation 

Y^ = st^l^+ f st^u{f{Y:)-dw:\ 

Jo 

understood in the classical Riemann-Lebesgue sense. Note that can be equivalently de- 
scribed as follows: with the expansion ([3]) of in mind, i.e. W = W(/3'), we have that 
VF" = W(/3"''), where, for each k > 1, (3"'''^ stands for the linear interpolation of (3^ with 
mesh ^. Therefore, it suffices to check that the conditions (i) and (ii) in Proposition 14.51 are 
satisfied by /3" := which is a matter of elementary computations (it can be also seen as 
a particular case of the forthcoming Proposition I4.8p . 

Together with Theorem ll.il we retrieve the following almost sure approximation result: 

Proposition 4.7. Under the hypotheses of Theorem M.li let be the Wong-Zakai approxima- 
tion of with mesh ^ and initial condition ip. Then, as n ^ cc, one has Af[Y—Y"'; C^{B^)] — t- 
a.s. 

This almost sure result in a non-linear situation is closely related to those of [5] or [2] , where 
Wong-Zakai' approximations for some parabolic type equations have been considered. We also 
note that convergence in law for this type of approximations in the framework of stochastic 
evolution equations has been studied in pU] . 

Now, let us turn to the proof of the weak approximation results of Theorem 11.21 and which 
successively involve the Donsker approximation /?"■ = and the Kac-Stroock approximation 
_ Qn^ both cases, we wish to exploit the criterion of Proposition 14.51 which naturally 
leads us to the following 2-step procedure: 

Step 1: Show that Condition (i) is satisfied, i.e., sup„E[|/3f - /S^pP] <Cp\t- sf. 

Step 2: Find a probability space {Cl,^,P), a sequence /3" and a Brownian motion (3, both 
defined on such that /3" ~ /3" and sup^g^,!] -^[1/3? - ^iP^] < CpU-^P for some fixed 

parameter v > 0. 

Once these two conditions have been checked, the proof of the weak convergence Y^ — )• Y 
in C^{Bry) becomes a straightforward consequence of Theorem 11.11 and Proposition 14.51 since 
W(;S"'') ~ W(/3"'') and accordingly, if Y^ denotes the solution of ([6]) associated with := 
W(^"'-), it holds that ~ y". 

Note that for both approximations S"" and 6"', the result in Step 2 will be derived from 
a Skorokhod embedding argument (see [27]). In the Donsker situation (Proposition 14. 9( ). this 
relies on a classical strategy towards the celebrated invariance principles (see \22\ Section 5.3]). 
In the Kac-Stroock situation (Proposition 14. IT]) . we will take advantage of an identification 
result due to Griego, Heath and Ruiz-Moncayo (see |14j). 

4.3. Donsker approximation. Here, we proceed to tackle the above 2-step procedure for 
the Donsker approximation S"". 

Step 1 (Donsker case): 



I ^ -h i 
n ' n 
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Proposition 4.8. For every p > \, there exists a positive constant Cp such that, for all 
< s < t < 1, 



snp E[\S^ - S^M <Cp\t 



Proof. First, note that 5" can also be expressed as 

ft 



(59) 



1=1 ''^ 

Then, by Lemma l4.H we have 

i=i -^^ 
< CpnPE[\Zi\^P] 

ft ^ ft 



u)duY < Cp\t-s\P. 



□ 



Step 2 (Donsker case): 



Proposition 4.9. Let (Z,)j(=N be a sequence of i.i.d. centered random variables with unit 
variance. Then, there exists a probability space {(},J-,P), a Brownian motion (3 defined on it 

and, for each n > 1, a family of independent random variables (^i"^)i=i,...,n with the same law 
as Zi, such that the following is satisfied. Set 

j ^ 



5r 



n 



+ 



1/n 



t G 



n n 



, with i G {1, . . . , n}. 



Then, for every integer p > 1, 



sup E[\^{t) - S^l'^P] < Cpn-P/\ 
te[o,i] 

Proof. As mentioned earher, it is based on a general Skorokhod embedding theorem (see [271 
p. 163]), which, in our particular situation, can be stated as follows : there exists a probability 

space (0, T,P), a Brownian motion ^ defined on it and, for each n G N, a sequence {Tj^"^}i=i,...,n 
of independent and positive random variables such that the random vector 

(^(rj")) , /3(rJ") + )) , . . . , ^(rj") + • • • + t^")^ 

has the same law as 

(^Zi,-^{Zi + Z2), . . . , ^(Zi + • • • + Z„)). 



Moreover, it holds that .^[rf"''] = E\{Zil^/nf\ = ^ and 



E\\r 



Set Tq := and T,j^^'^ := Yl]=i '^j^^' ^'^^ i > 1. With this notation, we can infer that 



in) I 



< CmE 



An) 



I \ 2m 



n 



< 



Crr 



for any m G N. 



n 
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We define now 
and 



i-l 



V- 7(n) t-{i-l)/n -(„)-. 

i=i 



if t G 



i — 1 z 



n n 



Observe that, if t G 



j— 1 i 



and hence, for t G 



i—l i_ 
n ' n 



we have that 



E 



\m - s\ 



< CpE 



< CpE 



1/n 

\m-~P{Tt\)\'' 
\m-i3{Tn)\'' 



+ CpE 



+ Cpn~P. 



,2p 



n 



Thus, we only need to bound the first term in the latter expression, and to this end, we will 
use the following decomposition: 



E 



^2 ' 



with 
and 



{|t-7;'"ll<n-l/n 



A^ = E 



On the one hand, the maximal inequality for Brownian motion yields 

+ E 



A'l<E max \P{s) - ^{t)\^P 

Lse[(t-n-i/4)V0,i] 



max |^(s) 

se[t, (i+n-l/4)Al] 



< 2E 



max <C„E\\^{n-'^/^)M <Cpn-P/^. 

L/iG[0,n-i/4] J 



On the other hand, by Cauchy-Schwarz inequality, we have 

^^<C,{£;[/3(i;(!i)^^]+i?[^(t)^^]}'^'{p(|i-Ta|>n~V4)} 

Note that, by Lemma |4.H 



1/2 



(60) 



(61) 



E 



— {Zi + ... + Z, 

n 



4p 



< CpU-'^Pi'^P < Cp. 



Thus, in order to estimate the term A2, we only need to study the probability appearing in 
([6T]) . To do so, observe first that since t G we have, for n such that | n~^^^ > ^ (that 

is, for n > 3), 

(62) 



i-l, 1 



Then, using again Lemma l4.ll we get 



n 



2 



p{\Tn 



i - 1 , 1 



> 



n 



i-l 



Therefore, A2 <Cpn which concludes the proof. 



(63) 
□ 
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4.4. Kac-Stroock approximation. Along the same lines as in the Donsker case, we proceed 
now to analyze the Kac-Stroock approximations based on 9"^. 



Step 1 (Kac-Stroock case): 

Proposition 4.10. For every integer p > 1, there exists a positive constant Cp such that, for 
aUO< s <t<l, 



sup 



EiK-e'^fp] < Cp\t-s\p. 



(64) 



Proof. We have that 



CpuPE 



CpuPE 



2p 



2p 



2p 



'-{ni<W2<-<M2p 



2p 



where in the latter equality we have used the symmetry of the integrand. Taking into account 
that the two possible values of random variable (— l)^(™i)"' i-^(""2p) only depend on the fact 
that the exponent is even or odd, we can write the latter expression above as 



Cprf 



'-{■Ul<«2<-<«2p 



|^((-l)^-i^(""2«)-^(""2-i))dni • • • du2p 



(65) 



Using that for ui < n2 < • • • < U2p, the random variables N{nu^J — N {nu^^_-^) are independent 
with Poisson distribution of parameter n{u2i — U2i-i), we have that ([65]) is equal to 



CpTf 



^{•ill<«2<---<U2p}^ 



-2n 



[ELi(«2,-"2,-i) 



dui ■ ■ ■ du 



2p 



This term can be bounded by 



CpnP 



'-{ui<U2} ' ' ' -'-{M2ij-l<'"2p}^ 



-2n 



[ELi(«2>-'"2«-i' 



dui ■ ■ ■ du2p 



CpUP 



<Cp{t 



t fU2 



-2n(«2-ni)^^^ ^\ ^ (J p 



s J s 
P 



2n 



(1 



-2n(u2— s) 



)dU2 



This concludes the proof. 
Step 2 (Kac-Stroock case): 



□ 



Proposition 4.11. There exists a probability space (0,^, P), a Brownian motion /3 defined 
on it and, for each n E N, a process ^" with the same law as 0" in such that, for any 
V G (0, |) and any p G N, 



sup E 

te[o,i] 



|/3(t)-r(t)|2f 



-pv 



(66) 



for some constant C, 



p,i/ ■ 



Proof. First of all, it is clear that we can suppose p(l/4 — u) > 1 (otherwise, we can use 
Jensen's inequality). Then, following the lines of |141 Section 2], we consider a probability 
space {Cl,T,P) with the following mutually independent objects defined on it: 

(i) a Brownian motion /3, 
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(ii) for each n € N, a sequence of independent random variables i € N} such that 

^1"^ has an exponential distribution with parameter 2-y/n, 

(iii) a sequence {fcj, i G N} of independent random variables such that P{ki = 1} = P{ki = 
-1} = 1/2. 

The i.i.d. random variables /ci^j"^ ^2^2"^ • • • ' verify i?[A;j^|"^] = and i?[(A;j^|"^)"] = 
Therefore, with the same result of Skorokhod as the one quoted in the proof of Lemma [4. 91 (see 
[271 P- 163]), there exists a sequence of independent positive random variables t[^\ Tg"^ . . . , 

such that ^(rf ^(r{"^+rf ^) . . . have the same law as ki^^\ ki^^'^ +k2^^^ ■ ■ ■ , respectively. 
Moreover, it holds that 



and, for each m £ N, 



E[[ 



<CmE\kiCy 



Set r/"'* := Yl)=i ''"j"^ define 



Then, let 9^ = {9^(t), t > 0} be a piecewise linear process given on the grid ri^"\ r2^"\ . . . by 



and ^"(0) = 0. The t-^^^s are independent random variables exponentially distributed with 
parameter 2n, and it is proved in [H] that the process 9" thus defined has the same law as 9^. 

Now, to show (j66p . we decompose the term E \j3{t) — as the sum of the following two 

terms: 



and 



E? := E 



E^ := E 



(n) 



Let us first study E^. If t belongs to := [T^_|, ) for some £ = 1, . . . , 8n, we have that 



in) 



So 



8n 



E^ = Y^E[\m-0"it)\''liteAr^ 



i=i 



8n 



< C,Y^E[\P{Ti'l\)-m\''Mt<^A^}]+CpY.^[\^^^^^^^ 



8n 



i=l 
8n 



i=l 



(67) 
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where, for the last inequahty, we have used the fact that f3{T^'^^) — f3{T^^-^) ~ ki^i. Now, for 
any fixed i £ {I, . . . , 8n}, write 



+ E 



mTtl) - ^(*)l'''l{*eAnl{|i_T(:) |>n-i/4}J • (68) 



|2p 



The first term in ()68p can be bounded with the same argument as in the proof of Proposition 
9] (see ([60]) ). which gives 



E 



As far as the second term in ()68p is concerned, we have 



E 



< {E[\m - ^{T^ltp^^'^' 



P(t G^^",|t-rg| >n 



-1/4 



1/2 



(69) 



and since 



E 



< c, 



pi 



we only have to focus on the probability appearing in ()69p . To do so, let us notice that 

p(t^A^,,\t-Tt\\>n-^/') 

< P{\t-fi'l\\<r^,\t-T^l\\>n-'/^) 

< P{\t - < rr,rt^ < \n^l\ \t - T^l > n^l^) + P(?f ^ > ^n^l"^) 

< P{\t - < It - > n"i/4) + Cpn-3p/2 



Then 



/-I 



2p 



+ E 



\T, 



in) 



I - 1 

2n 



2p 



n 



-1 2n 

where we have used the same argument as in ()63p to get the last bound. Going back to ([& 
we deduce that E^ < Cpnn~P^^ < Cpn~'^P, since p is assumed to satisfy p(| — z^) > 1. 

Eventually, we must deal with • fact, we have that 



-p/2 



E^ < {E[\e-it) - m\''] y^\p{t > ft^)}'^' < Cp{p{t > 1^ 



1/2 



where we have used Lemma 14.101 If we denote by Nn a Poisson process with intensity 2n, we 
can write 



8n. 



8n 



i=i i=i 
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~(n) 

because the rj 's are independent random variables exponentially distributed with parameter 
2n. The latter probability can be bounded by using Stirling's inequality, as follows: 

P{NM) > 8n) = f e--i^ < Ce"- f 

k=8n k=Sn ^/2^^k( j 

This lets us conclude the proof. □ 
5. Appendix A: fractional Sobolev spaces 



We gather here some classical properties of the fractional Sobolev spaces [Ba;^ 
which are extensively used throughout the paper. We recall the notations Ba for Ba,2 and B 
for Bq. Let us first label the following well-known regularizing properties of the semigroup (see 



Proposition 5.1. Fix two parameters A < a G M. Then, for every ip £ Bx and t > 0, 

\\SM\B^<ct~^''^^^\\ip\\B, , \\ASM\B^<ct~'-^''-^^\MB,. (70) 
and for every 'ip £ Ba, 

ll^iV'-V'llB. <ct""^||V||5. , ||A5<V||e. <ct"i+(""^)||V||B„. (71) 

The next results are taken from the exhaustive book [26] on fractional Sobolev spaces. With 
the notations of the latter reference, our space Ba,p (a £ G N) corresponds to F^^. Let 
us first report some properties regarding pointwise multiplication of functions. Due to the 
multiplicative perturbation in it is indeed natural that these results should intervene at 
some point. In the statement, the notation E ■ F C G must be understood as \\<p ■ ipWc < 
c IIv'II-bII^IIf for every G E,ip £ F. 

Proposition 5.2. The following properties hold true: 

(1) Section 2.4.4]; One has 

U-{0,l)cB^a ^f «>f -7, (72) 

2r 4 

and in particular: 

LP{0,1) ■ B C B^a ^/ «> ^- 

(2) Section 4.6. 1]^ Let ai < 02 be such that ai + 02 > and a2 > \ - Then 

Ba, ■ Ba, C Ba,. (73) 

In particular, Ba is an algebra as soon as a > j. 

(3) (1261 Section 4.8.2]; Let a > and pi,p2,p > 2 be such that 2a < j- (i £ {1,2}) and 
J- + J_ = 1. Then 

Pl P2 V 

Ba,pi ' Ba,p2 C Ba,p- (74) 

Let us also label here the classical Sobolev embedding 

Ba,pCB^ if 2a>-, (75) 
P 
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which yields in particular: 



Ba C Boo as soon as a > -. 

4 



(76) 



Finally, in order to handle the non-linearity in (jH), we resort at some point to the follow- 
ing stability result for composition of functions (see [261 Section 5.3.6]): if / : M — >■ M is 
differentiable with bounded derivative, then for every a G [0, ^] and ip Ba, 

||/(</')||b. <cH1 + II<^IIb.}- (77) 
Here, / is also understood as its associated Nemytskii operator, i.e., f{ip){C) ■= fi^i^,))- 



6. Appendix B: A priori estimates on the solution 

It only remains to prove the two a priori controls (j42p and (j43p for the solution y of (jH) (or 
equivalently (|24p ). To do so, we will rely on the following result, taken from jll) Lemma 6.5], 
and which extends the classical Garsia-Rodemich-Rumsey in two directions: 1) it covers the 
case of 5-variations and 2) it applies to more general processes defined on the 2-dimensional 
simplex ^2 = {{t, s) G [0, T]'^ : t > s}. 

Lemma 6.1. Let 5* = 6 or 6. For every a,f3 >0 and p,q > 1, there exists a constant c such 
that, for any R : S2 ^ Ba,p, 



N[R-,C^^{[Q,T]-Ba,r>)]<c[Up^^^^^^^^^^^ 



where 



VU 1 1 Ba 



0<u<v<T 



\v — u[ 



dudv 



Proof of Lemma \3.1(A In both cases, we will resort to the previous Lemma, which essentially 
reduces the problem to moment estimates. Thus, the following Burkholder-Davis-Gundy type 
inequality (borrowed from [21 Lemma 7.7]) naturally comes into play: for every a > 0, one 
has, by setting Uq := Q^/'^{B), 

ft 



E 



' St-u{f{Yu)-dWu) 

I3a 



E 



\St-u{f{Yu) 
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\HS{Uo,Ba:) 



du 



(78) 



where the notation HS(Uq, Ba) refers to the space of Hilbert-Schmidt operators defined on Uq 
and taking values in Ba- Note also that the family (A^Cfc) defines an orthonormal basis of Uq 
and accordingly 

1/2 



\St-u{f{Yu) • *)\\HSiUo,Ba) = Afe||5i„„(/(yj • efc) 



(79) 



Now, to show that Y G C'^^{Boo), observe first that for every q > 1 and any small e > 0, 



<CgE 



E[\\m 

< C,{e[\\ J'^ St^u{f{Y^ 



||2g 
tfllBi 

3 



dWu 



|2g 



+ E[(^l\\St^u{P ■ f{Yu) ■ fiY^mts^^duY'] }. 



(80) 



The second summand in ()80p is trivially bounded by Cp |t — sp"^^* since 

\\St-u{P ■ f{Yu) ■ nY^Ms. <c\t- u\--^-' IIP • /(y„) • f'{Yu)\\B <c\t- u\--i- 
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As far as the first summand in (|8U|) is concerned, observe that 

V \k\\St-u{f{Yu) ■ ek)\\l, < c V Afc |t - u\--^-^' \\f{Yu) ■ ekh < c|t - m| 

k ^ k 

which, owing to ([78]) and ([79]) . entails that 

ft 



-1-2. 



E 



St^u{fiYu)-dWu)\\l\ 



Going back to ([80]) . we are in a position to apply Lemma [6?T] and assert that Y £ C* (Boo) C 
C^^iBoo) (we recall that rj is assumed to belong to (0, |)). Note that since G B^, these 
computations also prove that sup^grQ j-i -BrUlfUg^ ] < oo for e > small enough, which will 

be used in the sequel. 

In order to show that Y G C^{B'y), let us write, like in ([80|) . 



E 



< c,{e[\\ J' St^u{f{Yu)-dWu)\ 

< c,{e[\\ St^u(.f{Yu)-dWu)\ 



+ E 



+ 



f\st^u{P ■ f{Yu) ■ fiYumB^duY'] } 



\t — u\ du 



2q 



(81) 



Then one has successively 



E 



\St-u{f{Yu) • *) 11^5(1/0,67) 
-2g{7-»?) 



< c,\t- ur'^''^-'">E[[Y,^k\\f{Yu) ■ efclll^ 



< Cq\t — u\ 



-2(}{7-r)) 



E[{Y.^'^y^^-^\\B,Jek\\l^y] (use m 



< c,\t-ur^''^^~^'^{y^{Xk-k*'^)y{l+ sup E[\\Yt\\l\ ]} (nsem), 
^ te[o,T] 



with 2(7 - r/) < 1 and Y^^i^k ■ k^V < ^- Thanks to ([78]) and ([79]), we can go back to §3) 
and deduce that E[\\{6Y)ts\\l^ ] <Cg\t- s\^'^' for some small e > 0. Since ijj G B^, this proves 
in particular that Y G C^{B^). 

Let us now turn to and notice first that 6K^ = L^6f{Y), so 



\\i6K 



Jtus lis 



< 



LZ\\c(B,B)\m{Y))us\\B 



Besides, since 



< cwj\t-u\^~''\\{6Y)us\\B (use (1321)) 

< cwj \t - ul^'"^ ^\\{6Y)us\\b + llaus^slle} 

< cw,f{\t - s\-^^''M[Y;C^\Boo)] + \t- s|^"''+^AA[y;CO(S-,)]} . 



St-u{Sif{Y))^, ■ dWu) + / St-uiP ■ f{Yu) ■ f'{Yu)) du, 



l-^-tslle 



<Cq\t- (use (HE]) and ([79]) as above). We are 



it is easy to check that E 

thus in a position to apply Lemma l6. II and conclude that K^^ G C2^^{B). 



STRATONOVICH HEAT EQUATION 



31 



□ 

Acknowledgments 

We thank Samy Tindel for helpful discussions on the topic of the paper. We are also very 
grateful to two anonymous referees for their careful reading and suggestions, which have led 
to significant improvements in the presentation of our results. 

Maria Jolis and Lluis Quer-Sardanyons are supported by the grant MCI-FEDER Ref. 
MTM2009-08869. 

References 

[1] G. Bal. Convergence to SPDEs in Stratonovich form. Comm. Math. Phys., 292(2) :457-477, 2009. 

[2] V. Bally, A. Millet, and M. Sanz-Sole. Approximation and support theorem in Holder norm for parabolic 

stochastic partial differential equations. Ann. Probab., 23(l):178-222, 1995. 
[3] X. Bardina, M. Jolis, and L. Quer-Sardanyons. Weak convergence for the stochastic heat equation driven 

by Gaussian white noise. Electron. J. Probab., 15mo. 39, 1267-1295, 2010. 
[4] H. Brezis. Operateurs maximaux monotones et semi-groupes de contractions dans les espaces de Hilbert. 

North-Holland Publishing Co., Amsterdam, 1973. 
[5] Z. Brzezniak and F. Flandoli. Almost sure approximation of Wong-Zakai type for stochastic partial differ- 
ential equations. Stochastic Process. Appl., 55(2):329-358, 1995. 
[6] R. Buckdahn and J. Ma. Stochastic viscosity solutions for nonlinear stochastic partial differential equations. 

I. Stochastic Process. Appl, 93(2):181-204, 2001. 
[7] R. Buckdahn and J. Ma. Stochastic viscosity solutions for nonlinear stochastic partial differential equations. 

H. Stochastic Process. Appl, 93(2):205-228, 2001. 
[8] R.A. Carmona and J-P. Fouque. A diffusion approximation result for two parameter processes. Probab. 

Theory Related Fields, 98(3):277-298, 1994. 
[9] G. Da Prato and J. Zabczyk. Stochastic equations m infinite dimensions, volume 44 of Encyclopedia of 

Mathematics and its Applications. Cambridge University Press, Cambridge, 1992. 
[10] A. Deya. A discrete approach to Rough Parabolic Equations. Electron. J. Probab., 16:1489-1518, 2011. 
[11] A. Deya, M. Gubinelli, and S. Tindel. Non-linear rough heat equations. Probab. Theory Related Fields, 

153(l-2):97-147, 2012. 

[12] C. Florit and D. Nualart. Diffusion approximation for hyperbolic stochastic differential equations. Stochastic 

Process. Appl, 65(1):1-15, 1996. 
[13] P. Friz and N. Victoir. Multidimensional dimensional processes seen as rough paths. Cambridge University 

Press, 2010. 

[14] R.J. Griego, D. Heath, and A. Ruiz-Moncayo. Almost sure convergence of uniform transport processes to 

Brownian motion. Ann. Math. Statist., 42:1129-1131, 1971. 
[15] M. GubineUi. Controlling rough paths. J. Funct. Anal, 216(1):86-140, 2004. 
[16] M. GubineUi and S. Tindel. Rough evolution equations. Ann. Probab., 38(l):l-75, 2010. 
[17] Y. Hu and D. Nualart. Stochastic heat equation driven by fractional noise and local time. Probab. Theory 

Related Fields, 143(1-2) :285-328, 2009. 
[18] M. Kac. A stochastic model related to the telegrapher's equation. Rocky Mountain ,J. Math., 4:497-509, 

1974. Reprinting of an article published in 1956, Papers arising from a Conference on Stochastic Differential 

Equations (Univ. Alberta, Edmonton, Alta., 1972). 
[19] I. Karatzas and S.E. Shreve. Brownian motion and stochastic calculus, volume 113 of Graduate Texts m 

Mathematics. Springer- Verlag, New York, second edition, 1991. 
[20] R. Manthey. Weak convergence of solutions of the heat equation with Gaussian noise. Math. Nachr., 

123:157-168, 1985. 

[21] R. Manthey. Weak approximation of a nonlinear stochastic partial differential equation. In Random partial 
differential equations (Oberwolfach, 1989), volume 102 of Internat. Ser. Numer. Math., pages 139-148. 
Birkhauser, Basel, 1991. 

[22] P. Morters and Y. Peres. Brownian motion. Cambridge Series in Statistical and Probabilistic Mathematics. 

Cambridge University Press, Cambridge, 2010. 
[23] A. Pazy. Semigroups of linear operators and applications to partial differential equations, volume 44 of 

Applied Mathematical Sciences. Springer- Verlag, New York, 1983. 
[24] M. Pinsky. Differential equations with a small parameter and the central limit theorem for functions defined 

on a finite Markov chain. Z. Wahrschemlichkeitstheorie und Verw. Gebiete, 9:101-111, 1968. 



32 



AURELIEN DEYA, MARIA JOLIS AND LLUIS QUER-SARDANYONS 



[25] H.P. Rosenthal. On the subspaces of {p > 2) spanned by sequences of independent random variables. 

Israel J. Math., 8:273-303, 1970. 
[26] T. Runst and W. Sickel. Sobolev spaces of fractional order, Nemytskij operators, and nonlinear partial 

differential equations, volume 3 of de Gruyter Series m Nonlinear Analysis and Applications. Walter de 

Gruyter & Co., Berlin, 1996. 
[27] A. V. Skorokhod. Studies in the theory of random processes. Izdat. Kiev. Univ., Kiev, 1961. in Russian. 

English Translation: Addison- Wesley, 1965. 
[28] Hector J. Sussmann. On the gap between deterministic and stochastic ordinary differential equations. Ann. 

Probability, 6(1):19-41, 1978. 
[29] G. Tessitore and J. Zabczyk. Wong-Zakai approximations of stochastic evolution equations. J. Evol. Equ., 

6(4):621-655, 2006. 

[30] S. Tindel. Diffusion approximation for elliptic stochastic differential equations. In Stochastic analysis and 
related topics, V (Silivri, 1994), volume 38 of Progr. Probab., pages 255-268. Birkhauser Boston, Boston, 
MA, 1996. 

[31] S. Tindel. Stochastic parabolic equations with anticipative initial condition. Stochastics Stochastics Rep., 
62(l-2):l-20, 1997. 

[32] J.B. Walsh. A stochastic model of neural response. Adv. m Appl. Probab., 13(2):231-281, 1981. 

Aurelien Deya: Institut Elie Cartan Nancy, B.P. 239, 54506 Vandoeuvre-les-Nancy Cedex, France. Email: 
Aurelien . DeyaSiecn . u-nctncy . f r 

Maria Jolis and Lluis Quer-Sardanyons: Departament de Matematiques, Facultat de Ciencies, Edifici C, 
Universitat Autonoma de Barcelona, 08193 Bellaterra, Spain. Email: querOmat . uab . cat , mj olisOmat . uab . cat 



